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A single spin-i /2 particle obeys the Dirac equation in d > 1 spatial dimension and is bound by 
an attractive central monotone potential which vanishes at infinity (in one dimension the potential 
is even). This work refines the relativistic comparison theorems which were derived by Hall [I|]. 
The new theorems allow the graphs of the two comparison potentials 14 and 14 to crossover in 
a controlled way and still imply the spectral ordering Ea < Es for the eigenvalues at the bottom 
of each angular momentum subspace. More specifically in a simplest case we have: in dimension 
d = 1, if Jq (Vi,(t) — Va(t))dt > 0, X € [0, 00 ), then Ea < Ei,-, and in d > 1 dimensions, if 
fgiVbit) — Va{t))d^'‘‘‘^dt > 0, r e [0, 00 ), where kd = t (j + and r = ±1, then Ea < Eb- 
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I. INTRODUCTION 

The comparison theorem of quantum mechanics states that if two comparison potentials are ordered, i.e. 14 < 14, 
then the discrete energy eigenvalues are ordered as well Ea < E},. In the nonrelativistic case this is a straightforward 
consequence of the min-max variational characterization of the discrete part of the spectrum 1 , |^ . In the relativistic 
case the Hamiltonian is not bounded below, and a variational analysis is more complicated y-l6|. However, comparison 
theorems have been established by other means in d = 1 and d = 3 dimensions [l[, in d = 2 dimensions Q, and in 
d dimensions Q, most recently by monotonicity arguments [ol-fT^. In Ref. [l^ Semay used the Hellmann-Feynman 
theorem [l^ to established a general comparison theorem for the Schrodinger and Dirac equations. 

In this paper we derive refined comparison theorems which allow the graphs of the comparison potentials to cross 
over in a controlled fashion and still imply definite ordering of the respective eigenvalues at the bottom of each 
angular-momentum subspace. This idea was first explored by Hall et al for nonrelativistic problems in d = 1 and 
d = 3 dimensions [l5l| and in d > 1 dimensions [l^ . and applied to Sturm-Liouville problems in El- In the simplest 
case one derives the spectral ordering Ea < Eb from the weaker potential assumption Ua < Ub, where Ui = Vi(t)dt, 
i = a or b. Since these refined nonrelativistic results were obtained without the use of a variational characterization 
of the discrete spectrum, similar reasoning could be applied to derive a basic relativistic comparison theorem for the 
Dirac equation [l|. The principal aim of the present paper is to go further and derive refined comparison theorems 
also for the Dirac spectral problem itself. 

In dimension d = 1, the energies compared are simply the lowest discrete eigenvalues. In d > 1 dimensions, the 
energies are the lowest eigenvalues in each angular-momentum sector. The derivations rely on a priori knowledge of 
the nodal structure characterized for central fields in Refs, [isl - l^ . We found it necessary to discuss the cases d = 1 
and d > 1 separately and to treat a small number of distinct classes of attractive monotone potentials. Sharper energy 
bounds can be obtained if the component wave functions are also known for the chosen base comparison potential. 
Simple sufficient conditions are derived in corollaries to the comparison theorems to make their use more immediate 
and straightforward. The results are illustrated by some specific examples. 
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FIG. 1: Ground state wave functions and (^ 2 , corresponding to the exponential potential V = —pe with P = 0.9, 
b = 0.5, m = 1, and energy is i? = 0.49233. 


II. DIRAC EQUATION IN ONE DIMENSION 

The Dirac equation in one spatial dimension for a single spin- 1/2 particle of mass m in natural units h = c = \ may 
be written [2l|: 


ai — - {E - U)cr3 + mj Ip = Q, 

where cti and 173 are Pauli matrices and the discrete energy eigenvalue E such that —m < E < m, The vector 

potential V (the time component of a four-vector) satisfies 

(i) V is even, i.e. V{x) = U(—x); 

(ii) V is nonpositive and bounded, i.e. Vo <V <Q, where Vb = P(0); 

(Hi) V vanishes at infinity, thus lim V = 0; 

fC—)-±oo 

(iv) V is attractive, that is monotone nondecreasing on [ 0 ,oo). 

By taking the two-component Dirac spinor as ip ) the above matrix equation can be decomposed into a 


system of first-order linear differential equations [mil: 


(fi'i = -{E + m- V)(fi 2 , 
^ 2 = iE-m-V)(pi, 


(la) 

(lb) 


where prime ' denotes the derivative with respect to x. For bound states, (fi and (^2 satisfy the normalization condition 


00 

(‘/5l, V^l) + {‘P2,¥>2) = j {^1+ ‘■pDdx = 1. 


For the reason which will be clear later, the comparison theorem below requires knowledge concerning the ground 
state. According to the Nodal Theorem of Ref. HI, the upper and lower components of the Dirac spinor, tpi and ip 2 
respectively, have definite and opposite parities and n 2 = ni + 1, where n^, * = 1 or 2, the corresponding number of 
nodes of ifi. Thus in the state with the smallest number of nodes, the upper component tpi is even and the lower one 
ip 2 is odd. From now on, without loss of generality, we consider the interval [0, 00 ) and assume that both components 
of the Dirac spinor lie above the x-axis, i.e. > 0 and (^2 > 0 on [0,oo). Then it follows from (lla|) - (llbl) that 
on [0, 00 ), < 0 and :/72 ^ 0 near the origin then (/?2 < 0 at infinity. As an illustration we plot {ipi,ip 2 } for the 

exponential potential V = — [1^ (see Figure 1). 
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III. REFINED COMPARISON THEOREMS FOR THE DIRAC EQUATION IN ONE DIMENSION 


We compare two problems with symmetric potentials 14 and Vf, and ground state energies Ea and Eh for which the 
system (fTa1)- (fTbl) becomes respectively 


and 


Let us consider the combination of equations: 

((2Ql) (/?9f. - 

which after some simplifications becomes 

{^la>P2b)' - {‘P2a‘Plb)' = {'fla>Plb + V2a‘P2b) [{Eb - Ea) - (14 - 14 )] • 

Integrating the left side of the above expression by parts from 0 to oo and using the boundary conditions, (^io(O) = 
‘^ib(O) = 0 and lim ip 2 a = Ihu ^ 2 b = 0, we find /„ [{^pia^ 2 by — {'^ 2 a^ib)'] dx = 0. Then we integrate the right side 

x—¥oo ^ 

to obtain 


-{Ea+m - Va)ip 2 a, 

(2a) 

{Ea-m - Va)ipia, 

(2b) 

-{Eb + m- Vb)ip2b, 

(3a) 

{Eb-m - Vb)ipib- 

(3b) 

Tib - ^^T 2 a + dMllT’la, 

(4) 


nco poo 

(Eh-Ea) +‘P2a^2b)dx = {^la^lb + ^2a^2b)iVb - Va)dx. 

Jo Jo 


(5) 


It follows from the last expression that if the wave functions have no nodes, so that the integrands have constant 
signs, and the potentials are ordered i.e. 14 < 14, then Ea < Eh- This is the coiMarison theorem which was first 
proved in d = 3 dimensions in [ij, in d = 2 dimensions in Q, and in d dimensions in Q. Later, using the monotonicity 
concept, the comparison theorem was proved for higher dimension cases d > 1 for the Dirac equation in 0 and the 
Klein-Gordon equation in [I^ for all the excited states. 

The Dirac equation admits exact analytical solutions for very few potentials. The above theorem allows us to obtain 
upper or lower bounds for any eigenvalue with the aid of suitable comparison potentials. But the comparison potentials 
can not cross each other, because in that case the integrands of (O change sign. Similarly to the nonrelativistic case 
we now derive refined relativistic comparison theorems which allow the graphs of the potentials to crossover in a 
controlled manner so that spectral ordering is predicted. 

Theorem 1: The potential V satisfies (i)-{iii), —14 < 2m, and has area. Then if 

g{x) = [ (14(i) - Va{t))dt >0, X € [0, oo), (6) 

Jo 

we have Ea < Eh- 

Proof: We integrate the right side of dU by parts to obtain 

poo poo 

/ {TlaTlb + T 2 a‘P 2 b){Vb - Va)dx = {(filaipib + T 2 aT 2 b)g\’^ - / gi^laTlb + T 2 aT 2 bydx, 

Jo Jo 

where g{x) is defined by ([6]). Since g(0) = 0 and lim (fia = lim (p 2 a = 0, relation ([5]) becomes 

x—¥oo x—foo 

pOO pOO 

{Eb - Ea) / {ipiaTlb + T 2 aT 2 b)dx = - g{ipiaTlb + ip 2 aT 2 bydx. (7) 

Jo Jo 

In order to find the sign of {(fiiaTu + T 2 aT 2 by we write 

{TlaTlb + T2aT2by = Wiipibip2a + W2ipiaT2b, 


where 


Wi = -Ea + Eb - 2m + Va - Vb and W 2 = Ea - Eb - 2m - Va + 14- 
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Let ip 2 represent either ^p 2 a or Lp 2 b- Suppose that ip 2 reaches its maximum at some point Xc] thus on [OjCCc] 

and (^2 — 0 oo [xcoo). It follows from (SB that Ea — m — Va > 0 on [0, so Wi < Eb — Vb — 3m. Since —Vb < 2m 
and —m < Eb < m, Wi < 0 on [0, Xc]- Then on [a;c,oo) equation (fTbl) implies Eb — m — Vb < 0, which leads to 
Wi < —Ea — m + Va, energy Ea satisfies —m<Ea < m so Wi < 0 on [xc, oo). Therefore Wi < 0 on [0, oo). Similarly 
it can be shown that W 2 < 0 on [0,oo). Thus {(pia^ib + ^ 2 a^ 2 b)' < 0, on [0,oo). Therefore if g{x) > 0 relation ([7]) 
ensures that Ea < Eb, which result completes the proof of the theorem. 

□ 

If we know the exact behaviour of the comparison potentials we can state simpler sufficient conditions for spectral 
ordering: 

Corollary 1: If the potentials cross over once, say at xi, Va < 14 for x G [0, xi], and 

POO 

g{oo) = {Vb - Va)dx > 0, 

Jo 

then Ea < Eb- If the potentials cross over twice, say at xi and X 2 , xi < X 2 , 14 < Vf, for x € [0, xi], and 

rX2 

5 (^ 2 )= / {Vb-Va)dx>0, 

Jo 


then Ea < Eb- 

We note that such application of Theorem 1 via the Corollary 1 can easily be extended. For example, consider the 
case of n intersections, n = 1, 2, 3, ..., and suppose again that Va < Vb for the first interval x G [0, xi\. Suppose 
now that the sequence |I4 — Va\dx, i = 1, 2, 3, ..., n, of absolute areas is nonincreasing (if n is odd then 

1 |I4 —14|da: > 114 — 14 Idee ), consequently it follows that fg(Vb(t) — Va(t))dt>0,xG[0, 00 ), and we conclude 

Ea ^ Eb- 

Remark: we now also consider theorems which take advantage of the known wave functions for one of the two 
comparison potentials. The general concept here is that we use these known wave functions for one of the eigenproblems 
along with an assumed relationship between the two potentials, and from these conditions we predict bounds on the 
eigenvalues of the second problem. In each such theorem we choose the base comparison potential to be 14 where in 
an application i may be chosen to be either i = a oi i = b; oi course, changing the base problem will also reverse 
the energy inequality from lower to upper bound, or vice versa. Now we state the second theorem (which allows the 
bottom of the potential to lie below —2m): 

Theorem 2: The potential V satisfies {i)-{iv) and has ipu and t(p 2 i-weighted areas, if 

ki{x) = f {Vb{t) - Va{t))(pu{t)dt > 0 and k 2 {x) = 

Jo 

where i is either a or b, then we have Ea < Eb. 

Proof: We prove the theorem for i = b; lor i = a, the proof is the same. We integrate the right side of ([S|) by parts 
to obtain 


{Vb{t) - Va{t))(p2i{t)tdt >0, X G [0, 00 ), 


( 8 ) 


{Eb - Ea) / + T2aT2b)dx 


7 1 ; 

klifla + k2 - 

00 |■°° 

kip'ia + k2 ( —] 

L X 

0 Jo 

\ X / _ 


where fci(x) and k 2 {x) are defined by dH]) for f = b. The expression \knpia + = 0, because A:i(0) = /c 2 ( 0 ) = 0 

and lim (pia = lim ip 2 a = 0. The function (pia is positive and decreasing, thus ip'ia < 0 on [0, 00 ). We know 

X^OO x—¥oo 

that i^ 2 a > 0 on [0,Xc] and ip 2 a < 0 on [xc,oo). It follows that {ip 2 alx) < 0 on [xc,oo). From the assumption 
{iv) and (l2b)) we conclude ip'fa < 0, that is to say y}2a is concave on [0,Xc], so lies below its tangents lines: thus 
T 2 a{x) > x:p 2 a{x), which implies {ip 2 alxy < 0 on [0,Xc]. Therefore {:p 2 alx)' < 0 on [0,oo). Finally, if ki{x) and k 2 {x) 
are both nonnegative it follows from the above expression that Ea < Eb. This completes the proof. 

□ 

The second theorem is sronger because the potential difference AV = 14 — 14 is multiplied by the decreasing factor 
ipii in fci(x), and by tip 2 i in k 2 {x), i = a, b, and this allows AV to be even larger than in Theorem I and still imply 
the spectral ordering Ea < Eb. Similarly to Corollary I, but now with the :/?ji-weighted areas, j = I or 2 and i = a 
or b, we can state the following sufficient condition for spectral ordering: 
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FIG. 2: The laser-dressed potential Va full line and the exponential potential 14 dotted line. 


Corollary 2: If the potentials cross over once, say at Xi, 14 < Vb for x € [0, xi], and 

POO POO 

fci(oo) = / (Vb — Va)(piidx > 0 and ^ 2 ( 00 ) = / (V& — Va)(p 2 ixdx >0, i = a or b, 

Jo Jo 

then Ea < Eb- If the potentials cross over twice, say at Xi and X 2 , Xi < X 2 , Va < Vb for x G [0, Xi], and 


rX2 rX2 

ki{x 2 ) = / {Vb — Va)y}iidx > Q and k 2 {x 2 ) = / {Vb — Va)ip 2 ixdx > 0, i = a or b, 
Jo Jo 


then Ea < Eb- 


Corollary 2 can be generalized as well for the case of n intersections: if 14 < 14 on a; G [0, a;i] and the sequences 
|14 — Va\y^ijdx and |14 “ Va\^P 2 jXdXi * = 1, 2, 3, ..., n and j = a or b, of absolute areas are nonincreasing 

(and, if n is odd, |14 - Va\ipijdx > /“ |14 - Va\(pijdx and \Vb - Va\(p 2 jxdx > /“ |14 - 141^52^2:^2:), then 

Io{^b{t) — Va{t))(pij{t)dt > 0 and jQ{Vb{t) — Va{t))(p 2 j{t)tdt > 0 on a; G [0, 00 ), so, according to Theorem 2, we 
conclude Ea < Eb- 


An example 


To demonstrate Theorem 1 we choose the laser-dressed potential 14 
the form 


and the exponential potential I 4 in 


14 . = - 


{x2 + ^ 2 ) 1/2 


and 


= _Rp-b\x\ 


14 = -Pe 


with a = 0.61362, a = 0.62, jd = 0.8, and b = 0.41; see Figure 2. Thus 14(0) = —0.989 and 14(0) = —0.8, and taking 
TO = 1, the condition —Vq < 2m is satisfied. The graphs of I 4 and 14 intersect at a;i = 0.94437 and X 2 = 4.13782. 
Then we calculate areas A and B 


A = [ {Vb-Va)dx = 0AU56 

Jo 


and 


B = 


rX2 

/ (14 - Vb)dx = 0.11455. 


Since A > B we have g > 0 therefore according to Corollary 1 we should have Ea < Eb, which we verify by calculating 
accurate numerical eigenvalues, i.e. Ea = 0.45657 < Eb = 0.52332. 
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IV. DIRAC EQUATION IN d > 1 DIMENSIONS 

For a central potential in d > 1 dimensions the Dirac equation can be written [s^ in natural units /i = c = I as 

where H=y asPs + mj5-\-V, 
at 

S — 1 

where m is the mass of the particle, V is an attractive spherically-symmetric potential, which will be defined later, 
and {Os} and /3 are Dirac matrices, which satisfy anti-commutation relations; the identity matrix is implied after 
the potential V. For stationary states, some algebraic calculations in a suitable basis, the details of which may be 
found in Refs, (dll - l^ . lead to a pair of first-order linear differential equations in two radial wave functions {'^i, 
namely 


I ^[ = {jn + E- V)'tp 2 - -yi’i, (9a) 

I = {m-E+ V)'tlJi +—'>p2, (9b) 

V r 

where r = ||r||, prime ' denotes the derivative with respect to r, kd = t (j + ^^), r = ±1, and j = 1/2, 3/2, 5/2, _ 

We note that the variable r is sometimes written oj, as, for example in the book by Messiah [s^, and the radial functions 

are often written tpi = G and tp 2 = F, as in the book by Greiner [s^- For d > I, these functions vanish at r = 0, and, 

for bound states, they may be normalized by the relation 


(ipuipi) + (V’2,V'2) = J (i’l +i^l)dr = 1. 

0 

We use inner products without the radial measure r*- •' because the factor r 2 is already built in to each radial 

function. We shall assume that the potential V is such that there is a discrete energy eigenvalue E and that equations 
(|^ - (|9b| are the eigenequations for the corresponding radial eigenstates. Throughout this paper we will consider 
only potentials which vanish at infinity, thus the above system at infinity becomes 

r = {m + E)'>p 2 , ( 10 a) 

\ ip2 = ~ (10b) 

Let us assume that 'ipi > 0 before vainshing, then it follows from (llOal) that {m + E)'tp 2 < 0. Thus either m + E > 0 
and ip2<0 01 m + E<0 and ip 2 > 0. By considering equation (Il0b|) . the first case leads to '02 ^ 0 and m — E > 0, 
so —m < E < m. The second case leads to '02 < 0 and m — E < but this is the contradiction: it follows from 
m + E < Q that E < 0 and from m — E < 0 that E > 0. Therefore we conclude that if the potential V vanishes at 
infinity then the discrete energy E is such that —m < E < m. 


V. REFINED COMPARISON THEOREMS FOR THE DIRAC EQUATION IN d > 1 DIMENSIONS 

As in one-dimensional case we need to know some characteristics of the nodeless state of the Dirac coupled equations 
(l9a]l - (I9bp . It follows from the Nodal Theorem of Ref. [1^ that in the state with no nodes kd < 0, and either '0i > 0 
and 02 < 0 or 01 < 0 and 02 > 0 for r G [0, c»); so from now on without loss of generality we suppose kd < 0 and 
01 > 0 and 02 < 0 for r G [0, 00 ). In Figure 3 we present an illustration of a node free state. Using (I9a1l - (l9bll and 
following the same argument as in one-dimensional case, we can obtain the corresponding relation for two comparison 
potentials 10 and I 4 


{Eb-Ea) / (01a01b + 02a02h)dr = / (10 - 10) (01a01h + 02a026)d?’- (H) 

Jo Jo 

From equation dm) we can eventually recover the basic comparison theorem [ij to the effect that if the radial 
components of the Dirac spinor are node free and 10 < 14, then Ea < Eb- 
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FIG. 3: The Dirac radial wave functions ipi and ^2 at the bottom of an angular-momentum subspace labelled by j, corresponding 


to the Woods-Sajcon potential in the form V = — 
and the energy eigenvalue is E = 0.62317. 


- ;:r^, with u = 4, i? = 2, o = 1.2, r = —1, d — S, j = 3/2, m — I, 

1 + 


A. Bounded Potentials 


Here we suppose that the potential V in the Dirac equations (Eili-dlbl is such that 

(i) V is nonpositive and bounded, i.e. Vb < H < 0, where Vq = H(0); 

(a) V vanishes at infinity, thus lim V = 0; 

r—^co 

(in) V is attractive, that is monotone nondecreasing on [0,oo); thus, when V is differentiable, V' > 0. 


We shall first prove the following lemma, which characterizes the behaviour of the Dirac radial wave functions at 

the bottom of the angular-momentum subspace labelled by j: for these nodefree states [l^, kd = — 

where j = 1/2, 3/2, 5/2, ... and d = 2, 3, 4, _ 

Lemma 1: the Dirac radial spinor components ipi and '02 o-t the bottom of an angular-momentum subspace labelled 
by j, which satisfy for the bounded potential V are such that 



01 

rl^dl 


< 0 and 


02 




>0, r G [0, 00 ). 


Proof: Near the origin the system (ISiI)-(l9b| may be rewritten as: 

f 0 " = 01 , ( 12 a) 

[ 0^' = 02 +m^-{E- Vof^ . (12b) 

Solutions of these equations involve Bessel functions. Hence for small r we can approximate them by simple powers 
in the following form 

r 01 = cir'?0 (13a) 

I 02 = C2r'J0 (13b) 

where ci and C 2 are constants of integration and parameters qi and q 2 are positive since both wave functions must 
vanish at the origin. After substituting (I13al) - (ll3bl) into (l^ - (l9bl) and dividing one equation by the other we obtain 
the following relation 

{qi + kd){q 2 - kd) = {m + E - Vo){m - E + Vo)r'^, 
which in the limit as r approaches 0 reduces to 


{qi + kd)[q 2 - kd) = 0. 















Since kd < 0 and qi > 0, it follows from the above expression that qi = —kd- Then equation (I9bl) becomes: 

C2(g2 - kd)r'^^~^ = ci{m - E + Vo)r~^'^. 

Equating the powers of r we obtain 92 = 1 ~ Also one finds 

Cl _ 1 - 2kd 

C 2 m — E + Vq 

According to (ii) and {in), lim [in — E + V) = m — E > 0, meanwhile m — E + Vo < 0, thus the quantity m — E + V 

r-^ nr> 

changes sign exactly once for r > 0. The ratio ci/c2 < 0, which means that and 'ip 2 have opposite signs (this is 
in agreement with our assumption for the nodeless state, that is '01 > 0 and 02 < 0 on [0, 00)). Finally, we conclude 
that near the origin radial wave functions behave as 


r 0l=Clr-'=^ (14a) 

1 02 = C2r-'=‘^+b (14b) 

Now let us make the following substitution 0i = and 02 = then the system of equations (l^ - (l9bl) 

becomes 

( R[ = {m + E-V)rR 2 , (15a) 

I = E + V)— + ^^^^R2. (15b) 

1 r r 

According to (0 and (I15a|) . R[ < 0 which is equivalent to the lemma’s first inequality. Since m — E + V has to change 
sign from negative to positive, thus for large r, according to (|15bl) . i?2 > 0. In order to determine the behaviour 
of the i?2 near the origin, we expand Ri and i?2 in power series, i. e. i?i = oq + air + a 2 r^ + + ... and 

R 2 = bo + bir + b 2 r‘^ + + ..., then system (|15al) - (ll5bl) implies 


f i?i = tto + a2r^ + ..., (16a) 

1 i?2 = &o + + ..., (16b) 

where if oq > 0, then bo < 0, 02 < 0, and 62 > 0. Thus i?2 > 0 near zero. 

Let us suppose that the function R 2 is decreasing on some interval (ri, r2), i. e. i?2 < 0 on (ri, r2) and 
-^2(^1) = -^2(^2) = 0. Then it follows from the Rohe’s Theorem that there is at least one number Tc G (ri, r 2 ) such 
that R'-iirc) = 0, which is equivalent to 

({m — E + V)—-\ - - -i?2^ =0 at r = rc 

\ r r J 

or, using (Il5bl) . 

V'Ri + {m-E + V)R[ + {2kd - 2)R'^ =0 at r = r^. 

In the above expression first two terms are nonnegative and {2kd — 2 )R '2 is strictly positive, which yields a contradiction. 
Hence i?2 — 0 ^ ^ [Oj o®) corresponds to the lemma’s second inequality. 

□ 

Theorem 3: The potential V, satisfies {i)-{ii), Vq > —2m, and has -weighted area, if 

rj{r) = f iVb{t) r € [0, cxd), (17) 

^0 


then we have Ea < Et,. 

Proof: Let us integrate the right side of (ITT]) by parts in the following way 


I 

(14 - Va)('01a0’lb + i’2a^f2b)dr = (01a'01b + '02a02b)-5|^ 


'f’la'f’lb + 02a'02b y 

J 


dr, 
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where ri{r) is defined by (1171) . Since 77 ( 0 ) = 0 and lim tpi = lim '!/'2 = 0 with the respective asymptotic forms 

_ _ r—>-oo r—>-oo 

(|14al) - (ll4bl) . relation (ITTl) becomes 

(Eb-Ea) j {ipia'<pib + ^2ail^2b)dr = - J 7 ? ^^ dr. (18) 

Using (Pal) - (l9bl) we find 


/ V^laV'lb + •02a’02b 


/ 


= VFi 


l/’lbl/’2a 

I 


+ W2 


'01a'02fe 

J’^l^d I 


+ ^kd 


•02a'02fe 
1 + 1 


where 


lUi — 2m + Ea — Eb — Uj + Vf, and IU 2 — 2m — Ea + Eb E Va — Vb- 


We know that quantity m — E + V changes sign from negative to positive. When m — E + V <0we have Wi > 
3m — Eb + Vb- Since —m < E < m and V > —2m, then Wi > 3m — Eb + Vb > 0. When m — E + V >0itis 
straightforward that Wi > m + Ea — Va >0. Therefore Wi > 0 for r G [0, 00 ). Similarly it can be shown that 

"V'laV'lb + ■4’2a'tp2b\' 


W 2 > 0. Thus the derivative 
m we have Ea < Eb- 


n2|fcdl 


< 0 and according to the theorem’s assumption (inD and equation 


□ 


We note that above theorem, as well as Theorem 1, does not require a nondecreasing potential V on [0, 00 ), i. e. V 
can decrease on some intervals. As in the one-dimensional case, if we know more precise behaviour of the comparison 
potentials, we can state simpler sufficient conditions: 

Corollary 3: If the potentials cross over once, say at ri, and Va < Vb for r G [0, ri], and 

nOO 

77 ( 00 )= / {Vb-Va)r‘^^^'^^dr>0, 

Jo 

then Ea < Eb- If the potentials cross over twice, say at ri and r 2 , ri < r 2 , Va <Vb for r G [0, + 1 ], and 

^(^ 2 )= [ {Vb-ya)r^^^'^^dr>Q, 

Jo 

then Ea < Eb- 

We can extend the above corollary in the following way: assume that comparison potentials have n intersections, 
77 , = 1, 2, 3, ..., and 14 < 14 onr G [0, ri]. Also assume that |14— 14 | 7 '^l^‘^ldr, i = 1, 2, 3, ..., n is nonincreasing 

sequence (if ?7 is odd then JJ'"' J14 —14|r^l^‘^ldr > 114 — 14 I dr ), hence J^(yb(t) — Va{t))t'^^’^'^^dt>0,rG [0, 00 ), 

and we conclude Ea < Eb. 

Theorem 4: The potential V satisfies {i)-{iii), and has and if-weighted areas, if 

= [ (Vb(t) - 14 (t)) 7 />ii(t)tl'''*ldt > 0 and A 2 (r) = f (14(t) - 14 (t))(- 4 ’ 2 z(t))t''''^'^^dt > 0, (19) 

Jo Jo 

r G [0, 00 ), where i is either a or b, then we have Ea < Eb. 

Proof: We prove the theorem for i = b; for i = a, the proof is the same. After integrating the right side of (llip by 
parts we obtain 

pOO 

/ (14 - Va)(V’lalAlb + 'll’2a'll’2b)dr = 


\ tpla , . i’2a 
''l-iTTr + -^2- 




' 7*1 1 + 1 


Ai 


'f’la 

r4dl 


+ A 2 


'f’2a 
7 -l^dl + l 


dr, 


where Ai(r) and A 2 (r) are defined by m for i = b. Since Ai(0) = A 2 ( 0 ) = 0 and lim ifia = lim if 2 a = 0 with the 

_ _ r—¥oo r—>00 

respective asymptotic forms (ll4aD - (ll4bll . expression (fTT|l becomes 


poo 

{Eb - Ea) / {'Ipla'iplb + i’2ai’2b)dr = - 

Jo 


Ai 


r’l^d I 


+ A 2 


1p2a 
1 + 1 


dr. 
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FIG. 4: The graph of the integrand I 


Then Lemma 1 and the theorem’s assumptions ensure that Ea < Eh. 

□ 

Corollary 4: If the potentials cross over once, say at ri, Va < 14 for r € [0, ri], and 

poo pOO 

Ai(oo)= / {Vb-Va)'il^ur^^^^dr>0 and A 2 (oo) = / {Vb - > 0, i = a or h, 

Jo Jo 

then Ea < Eh- If the potentials cross over twice, say at ri and r 2 , ri < r 2 , Va < 14 for r € [0, ri], and 

Ai(r 2 ) = [ {Vb-Va)'ipur^^‘‘^dr>0 and \ 2 {r 2 ) = [ (14 - > 0 , i = a or b, 

Jo Jo 

then Ea < Eh- 

As before we can generalize Corollary 4 to allow n intersections, i.e. if 14 < 14 on r G [0, ri] and both sequences 
of absolute areas |(14 — 14 ) 4 ’ii?’'^‘^'|dr and |(14 — i = a or b, are nonincreasing (if n 

is odd then we assume |(14 - Va)ipiir^^‘‘^\dr > |(14 - 14)V’ii»''^'^'and |(14 - Va)'(p 2 ir^'^‘^^~''^\dr > 

fZ \irb-Va)'ilJ 2 ir^''‘‘^^^\dr), then integrals /J’(14(t)-14(t))'*/'ii(i)i'''''''^^ > 0 and /Q’'(lb(t)- 14 (t))(-V’ 2 i(t))t'''‘^'+^dt > 0 
for r G [ 0 , oo), thus Ea < Eh- 


An example 


As an example we consider Theorem 3, in particular Corollary 3 for the case of many intersections. We take the 
following comparison potentials 


14 = - 


ur^ + a 


1 + 


V sin(Kr^ + s) 


Atr"’ 


and 14 = — 


wr^ + b 


lia = P,a = s = b, u = K = w,T = —1, j = 1/2, and d = 3 the substitution z = ur^ + a transforms the integral (fTTl) 
into 


/ {Vh-Vaydt=Z^ 

lo 3m , 


smx 


dz. 


smz 


The integrand I = —— is plotted on Figure 4. Choosing a = 2.04, and calculating numerical values, we find that 


the first area is bigger then the second one: 

I sinzl 


dz = 0.0965 > 


sin z 


dz = 0.0962. 


.27r 
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The sin 2 is the periodic function, thus | sina;| 
it is clear that 



I sinyl where x £ [{k — l) 7 r, kir] and y 


sinz 


dz > 


^(fc+i)^ 


' kn 


Sin 2 


dz. 


a; + TT, fc = 3, 4, 5, ..then 


Therefore 

{Vb - Va)r‘^dt > 0 , 

because successive positive and negative areas of the integrand do not increase in absolute value. Thus rj > 0 and by 
Theorem 3 we have Ea < Eb- This prediction is verified by accurate numerical calculations: for a = 3.4, a = 2.04, 
V = 0.4, and u = 7 the comparison potentials intersect at infinitely many points (see Figure 5), with m = 1 so that 
condition —2m < Vq is satisfied. Accurate numerical eigenvalues are Ea = 0.99427 < Eb = 0.99542. 




FIG. 5: Left graph: potential Va dashed line and 14 full line. Right graph: function dashed line and full line. 


B. Unbounded Potentials 

f{r) 

We consider a class of unbounded potentials of the form V{r) = -, where the bounded factor / satisfies: 

r 

{iv) f is nonnegative and bounded, i.e. /o > / > 0 , where /o = lim /; 

r->-0+ 

(u) / is monotone nonincreasing on [ 0 , oo), so f' < 0 ; 

/ 

(vi) V vanishes at infinity, so lim — = 0. 

r —T 


V I I V 

For instance, for the Coulomb potential V = , the function f = v, for the Yukawa potential [39|| V = -t—, the 


function / = for the Flulthen potential V = -- 


,\r 




the function / = 


vr 


,\r 


- 1 


, and so on. 


Lemma 2: the Dirac radial spinor components "01 and 'ip 2 o-t the bottom of an angular-momentum subspace labelled 

f 


by j, which satisfy for the potential V = - are such that 

r 


V'l 


I / 


< 0 and 


i’2 

r\^d\ 


>0, r G [0, oc). 
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Proof: For small r analysis of the Dirac coupled equations yields the asymptotic forms: 




V'2 = '02 


j,2 

/O 


r 

r 


These are Cauchy-Euler equations with solution in the form [4l| : 

01 = cir^, 
02 = C 2 r'^, 


( 20 a) 

( 20 b) 


( 21 a) 

( 21 b) 


where Ci and C 2 are constants of integration, and the parameter 7 has to be positive because the wave functions must 
vanish at the origin. Substitution of (I21al) - (l21bl) into (l^ - (l9bl) yields 

7 = / 0 C 2 /C 1 - kd, 

7 = -foCilc2 + kd- 


The solution of the above system is: 7 = \/k‘d — Jq and — = 


C2 kd ± yjkl- 


. As in sec. A, 02 /ci < 0, which is in 


Cl fo 

agreement with our assumption for the nodeless state. Therefore near 0 the wave functions behave as 

f i/ji = , 

\ '02 = C2r'0^‘^~0 . 

We now substitute tpi = and 'ip 2 = r~^'^P 2 , into (|5a)) - (l9bl) to obtain: 

= (m + A - V)P2, 

, , , 2 fcd 

= (m - A + V)Pi + —P 2 . 


(24a) 

(24b) 


Clearly, P[ < 0 which is equivalent to the lemma’s third inequality. Near the origin P 2 behaves as thus 

P 2 > 0 near 0. At infinity (I24bl) becomes P 2 = {m — E)Pi, so P 2 > 0 near infinity. Let us assume that P 2 < 0 on some 
(ri, r 2 ), so P 2 (ri) = P 2 {r 2 ) = 0. Then by Rolle Theorem, there exists Vf. such that P^'irc) = 0, which corresponds to 


/ y\ p' 

(m- E + V)Pi+ iV+ -\ Pi + i2kd-l)y+ir 


p 

— P) — = 0 at r = 


/ . V f ,, , 

Since V = -, the expression V -\ -=-is nonnegative according to (v). Therefore in the above expression 

r r r 

the first two terms are nonnegative and the last two are strictly positive, which observation reveals a contradiction. 
Therefore P 2 > 0 on r G [0, 00 ) and this is equivlent to the lemma’s last inequality. 

□ 

Now we state and prove the refined comparison theorem for a special class of unbounded potentials: 

f 

Theorem 5: The potential V =- , where f satisfies {iv)-{vi), has and ~weighted areas, if 

r 

Pi(t)= [ {Vb{t)-Va{t))fiu{t)t^’^'‘^dt>0 and P 2 {r) = f (t4(t) - Va(t))(-02i(t))t'''‘^'dt > 0, (25) 

Jo Jo 

r € [0, 00 ), where i is either a or b, then we have Ea < Pb- 

Proof: We prove the theorem for i = a; for i = b, the proof is the same. As in sec. A we integrate the right side of 


m by parts to obtain 

^00 pQ 

/ (Vb - v;)(01a01b + i’ 2 a'>p 2 b)dr = - 
Jo Jo 


m 


01a 

^1 kd I 


772 


02a 

r\kd I 


dr, 
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where ? 7 i(?’) and 772 (r) are defined by (|25ll for i = a. Then it follows from last two inequalities of the Lemma 2 and 
(051) that Ea < Ef,. 

□ 


Corollary 5: If the potentials cross over once, say at ri, Va < 14 for r € [0, ri], and 

poo poo 

pi(oo) = / {Vb-Va)ipiir^'''^^dr >0 and ^ 2 ( 00 ) = / {Vb - Va){-ip 2 i)r^'''‘^dr > 0, i = a or b, 

Jo Jo 

then Ea < Eb- If the potentials cross over twice, say at ri and r 2 , ri < r 2 , Va < 14 for r € [0, ri], and 

Pi{r 2 ) = f (14 - > 0 and p 2 {r 2 ) = f (Vb - Va){-ip 2 i)r^’'‘‘^dr > 0, i = a or b, 

Jo Jo 

then Ea < Eb- 

The above Corollary can be generalized up to n intersections: if 14 < 14 on r e [0, ri] and both sequences of absolute 
areas |(14-14)4’ii’’'^‘^'M’’and |(14-14)4’2i7’'^‘^'M7', i = a OT b, are nonincreasing (if n is odd then we assume 
\iyb-Va)tjjur^''‘‘^\dr > /“ 1(14 - and |(14 - 14)7/'2*rl'=‘^l |dr > /“ |(14 - 14)7/'2jrl'=‘^l |dr), then 

Vilr) > 0 and 772 ( 7 ’) > 0 for r G [0, 00 ), thus Ea < Eb- 


An example 

Here we will demonstrate first part of Corollary 5, i. e. the case of one intersection. For the comparison potentials 
we choose the Hulthen potential 14 and the Coulomb potential 14: 

Q7 B 

14 = — - - T and 14 =-. 

gar _ I j. 

The above potentials intersect at exactly one point for a = 1/5, a = 3/10, and 4 = 0.508; Figure 6 , left graph. We 
can see that 14 < 14 before the intersection point. If pi ( 00 ) and ^ 2 ( 00 ) are both nonnegative, according to Corollary 
5, we have Ea < Eb. 

The solutions of the Dirac Coulomb problem are well known. In three dimensional ground state for j = 1/2, r = —1, 
and 777 = 2 the eigenvalue is Q 

Eb = 27 


and the wave functions are 

1 / 1 ,) /24(1±7)(W 

7/26 j Y r (27 + 1 ) e^hr ’ 

where 7 = -^/l — (3'^, and F is the gamma function (the wave functions are plotted on Figure 6 , right graph). Then 
pi(oo) = 0.00113 and ^ 2 ( 00 ) = 0.00031 and, according to accurate numerical calculation, Ea = 1.58604 < Eb = 
1.72271. 


C. Potentials less singular than Coulomb. 

We characterize this class of potentials in the following way: 

{vii) V is nonpositive and unbounded, i.e. H < 0 and lim V = — 00 ; 

--^0+ 

(via) V is attractive, that is monotone nondecreasing on [0, 00 ) so V' > 0; 

(ix) V vanishes at infinity, thus lim H = 0; 

r—foo 

(x) H is less singular than Coulomb potential, i.e. lim rV = 0. 

r—)-0+ 

V V V 

Examples of such potentials are: V = - , V = - j-, and V = - j—, 

'pQ pQ I pO pQ ^OT 

constants. Then, as before, we first prove a lemma, namely 


where q G ( 0 , 1 ), v and b are positive 
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FIG. 6: Left graph: The Hulthen potential Va dotted line and the Coulomb potential 14 full line. Right graph: Upper 
component of the Dirac spinor tpu dotted line and lower component 'i /)26 full line. 




FIG. 7: Left graph: Function Vatpur dotted line and Vbtpibr full line. Right graph: Function Va{—ip2b)r dotted line and 
Vb{—ip2b)r full line. 


Lemma 3: the Dirac radial spinor components ipi and rp 2 at the bottom of an angular-momentum subspace labelled 
by j, which satisfy for the potential V, which satisfies (i)-{iv), are such that 


f’l 

rl^dl 


< 0 and 


•02 y 

+ l J 


> 0 , 


r G [0, oo). 


Proof: Near the origin the above class of the potentials can be approximated hyV = 
after some rearrangements becomes: 


—, then the system (l^ - (l9bl) 


01 ( 

(kd{kd + 1 - g) 

v‘^ 

j,2 

j.2q 

02 ( 

f kd{kd - 1 + g) 


1 J.2 

j.2q 


/ q 




-02-- 

r 


(26a) 


(26b) 
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FIG. 8: Left graph: The graph of the integrand (Vi, — Va)'tl>ibr. Right graph: The graph of the integrand {Vb — Va){—'tl) 2 b)r. 


Solutions of these equations are given in terms of Bessel functions. Therefore for small r we can approximate them 
by simple powers 


ipi = Cir'Ji, 

V'2 = C2r'^", 


(27a) 

(27b) 


where ci and C 2 are constants of integration and parameters qi and 92 are positive since both wave functions must 

Cl 

vanish at the origin. Then, following the proof of Lemma 1, we find qi = —kd, 92 = ~ 9 + 1, and — = 

C2 

1 — q — 2kd 

-< 0. Thus, near the origin radial wave functions behave as 

V 

V'l = 

l/;2 = C2r-'=-^-«+^ 

Now let us make the following substitution i/ji = and ip 2 = then the system of equations (I9a1) - (j9bp 

becomes 


R'^ = {m + E — V)rR 2 , 

R 2 = {m — E + V) — -\—— -i? 2 . 

r r 


(29a) 

(29b) 


According to (vii) and (|29al) . R[ < 0 which is equivalent to the lemma’s first inequality. Near 0 the function i ?2 
behaves as C 2 r~‘^, thus R 2 > 0 near 0. Equations (I29al) - (l29bl) are exactly the same as (I15al) " (ll5bl) . thus it can be 
proved by contradiction that i ?2 > 0 on r G [0, 00 ) which is the same as the lemma’s second inequality. 

□ 

Since Lemma 1 and Lemma 3 have the same conclusions, the following theorem along with the corollary can be 
stated and proved as in the bounded case: 

Theorem 6: The potential V satisfies {vii)-(x), has and -weighted areas, if 

Ti{r) = [ {Vb{t)-Va{t))ipu{t)t^''''^dt >0 and fJ. 2 {r) = ( (V),(<) - > 0, (30) 

do do 

r G [ 0 , 00), where i is either a or h, then we have Ea < Eh- 

Corollary 6: If the potentials cross over once, say at ri, 14 < Vb for r G [0, ri], and 

poo poo 

fii{oo) = / {Vb-Va)'tpiir^’''‘^dr >0 and ^2(00)= / {Vb - Va){-'ip2i)r^'"'^^^^dr > 0 , 

Jo Jo 


i = a or 
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then Ea < Eb. If the potentials cross over twice, say at ri and r 2 , ri < r 2 , Va < 14 for r € [0, ri], and 

pr2 pr2 

Mi(^ 2 )= / {Vb-Va)'4)ur'''^'^^dr>0 and p. 2 {r 2 ) = / {Vb-Va){-'f 2 i)r''^'^''^^dr >0, i = a or b, 

Jo Jo 


then Ea < Eb- 

As before Corollary 6 can be generalized to the case of n intersections. 


D. General refined comparison theorem 


In this section we shall construct the refined comparison theorem which combines all the above cases, i.e. bounded, 
unbounded, and less singular than Coulomb potentials. First, we state the following lemma: 

Lemma 4: the Dirac radial spinor components ipi and rp 2 at the bottom of an angular-momentum subspace labelled 
by j, which satisfy for the bounded potential (satisfies (i)-(iii)), or for the unbounded potential (satisfies 

(iv)-(vi)), or for the potential less singular than Coulomb (satisfies (vii)-(x)) are such that 


i’l 

rl^dl 


< 0 and 


i’2 y 


> 0 , 


r G [0, oo). 


Proof: The proof for the bounded and less singular than Coulomb potentials follows from Lemma 1 and 3 respectively, 
so we have to prove the Lemma for the unbounded case only. 

Let us consider equations (Hi!)-® with substitution ifi = r ^‘‘Gi and 'ip 2 = r ^’‘~^^G 2 - 


G[ = {m + E - V)rG 2 , 

G '2 = {m-E + V)— + ^^±^G 2 . 

r r 


(31a) 

(31b) 


Then the lemma’s first inequality follows from (I31al) . We know from the proof of Lemma 2 that G 2 behaves near the 

origin as where the constant C 2 < 0, therefore G 2 > 0 near 0. Since {m — E + V) > 0 near infinity, 

then according to (j31bl) . function G 2 > 0 near infinity. Finally, it can be shown by contradiction that G 2 > 0 on 
r G [0, 00 ), which ends the proof. 

□ 

Using the above Lemma, the general refined comparison theorem can be proved: 

Theorem 7: The potential V satisfies either (i)-(iii), or (iv)-(vi), or (vii)-(x), has and if-weighted 

areas, if 


Ci{r)= [ iVb{t)-Vait))ijju{t)t^''''^dt>0 and C, 2 ir) = j (I4(t) - I4(t))(-V’2i(t))t''''^'^^dt > 0, 
Jo Jo 

r G [0, 00 ), where i is either a or b, then we have Ea < Eb- 

Corollary 7: If the potentials cross over once, say at ri, Va < 14 for r G [0, ri], and 

poo poo 

Ci(oo) = / iVb-Va)ipur^’''‘^dr>0 and ( 2 ( 00 )= / (Vb - Va){-if 2 i)r^'"''^~''^dr > 0, i = a or b, 

Jo Jo 

then Ea < Ei,. If the potentials cross over twice, say at ri and ^ 2 , ri < ^ 2 ; 14 < 14 for r G [0, ri], and 

pr2 pr2 

Ci(?' 2 ) = / (Vb — Va)ifiir^^'^^dr > 0 and C, 2 {r 2 ) = / {Vb-Va){-i> 2 i)r''^'^''^^dr >Q, i = a or b, 

Jo Jo 


(32) 


then Ea < Eb- 

As before Corollary 7 can be generalized to the case of n intersections. 
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FIG. 9: The Coulomb potential 14 dotted line and the Sech-squared potential 14 full line. 


An example 


In that section we will demonstrate the second part of Corollary 7, using unbonnded Coulomb potential 14 and 
bounded Sech-squared potential 14, which has been known under several different names since early days of quantum 
mechanics such as the Poschl-Teller potential [i^ or the Eckart potential (45 |: 

Va = —— and Vf, = —/3sech^br. 
r 


If a = 0.579, P = 0.3, and b = 0.2 then 14 and 14 intersect at exactly two points ri = 2.41742 and r 2 = 5.66301; see 
Figure 9. According to @, in d = 3 dimensions for j = l/2,r = —1, and m = 1 the wave functions and the eigenvalue 
are: 


^1“ \ = ± / 7) (2Qr)T' 

Tp 2 a J y r(27 + 1) e“’' 


and Ea = 7 , 


where 7 = \/l — and T is the gamma function. Then direct calculation shows that Ci(?' 2 ) = 0.18778 and C 2 (?’ 2 ) = 
0.00084. Thus Corollary 7 implies that Ea < Ei,, which agrees with the accurate numerical energy eigenvalies: 
Ea = 0.81533 <Eb = 0.88318. 


VI. CONCLUSION 

In this paper we have rederived the earlier comparison theorems in d dimensions and then refined these theorems 
by allowing the comparison potentials to crossover in a controlled manner. Because of the different form of the Dirac 
coupled equations in d = 1 and d > 1 dimensions we have studied these cases separately. We have shown that in one 
dimension the condition 14 < 14 , which leads to the spectral ordering Ea < Eb, can be replaced by weaker condition 
Ua < Ub, where Ui{x) = Vi(t)dt i = a or b. Since the potential cross-over conditions depend on the detailed 
behaviour of the radial wave-functions components, we have found it best in d > 1 dimensions to establish first a 
separate refined comparison theorem for each of a number of interesting classes of potential: Theorems 3 and 4 for 
bounded potentials. Theorem 5 for unbounded potentials, and Theorem 6 for the class of potentials less singular than 
Coulomb. Finally, we have summarized the refined comparison results in Theorem 7, which combines all of the above 
types of potentials and states that if Uia < Uib and C /20 < U 2 b then Ea < Eb^ where Uuir) = 
and U 2 i{r) = Jg Vi(t)(—'ip 2 j(i))t^^'^^~^^dt, where i, j = a or b. 
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FIG. 10: Left graph: The graph of the integrand (14 — Va)tpibr- Right graph: The graph of the integrand (14 — Va){—'<p 2 b)r^- 


For practical reasons we have also established weaker sufficient conditions, as corollaries, which guarantee in simple 
ways that the comparison potentials crossover so as to imply definite spectral ordering. These results are illustrated 
by some specific examples. 
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